High-frequency effects in the FitzHugh-Nagumo neuron model 
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The effect of a high-frequency signal on the FitzHugh-Nagumo excitable model is analyzed. We 
show that the firing rate is diminished as the ratio of the high-frequency amplitude to its frequency is 
increased. Moreover, it is demonstrated that the excitable character of the system, and consequently 
the firing activity, is suppressed for ratios above a given threshold value. In addition, we show that 
the vibrational resonance phenomenon turns up for sufficiently large noise strength values. 
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Non-linear noisy systems have been studied with ever 
growing interest due to its applicability to the modeUing 
of a great variety of phenomena of relevance in physics, 
chemistry and the life sciences Perhaps, the sim- 

plest non-hnear noisy system studied in the literature is 
the bistable system, which has been used successfully to 
illustrate the phenomenon of stochastic resonance. In 
addition, the spiking activity of a neuron has been the- 
oretically studied using nonlinear excitable noisy models 
y], among them, the FitzHugh-Nagumo (FHN) system 
|4| being one of the most utilized due to its simplicity. 
Besides, networks of FHN units have been considered as 
simplified models useful in the description of both the 
cardiac tissue 0, S 13 and reaction-diffusion chemical 
systems 

Recently, it has been shown, both theoretically and ex- 
perimentally, that the addition of a high-frequency (HF) 
signal results in an improvement of the stochastic res- 
onance in a bistable optical system PJ. Nevertheless, 
many experimental studies have suggested that HF (non- 
ionizing) fields may damage several structural and func- 
tional properties of neuronal membranes in single cells, 
as well as to cause a number of negative physiological 
effects on typically excitable media such as the cardiac 
tissue Besides, it has been found that certain HF 

signals are able to suppress the steady directed motion 
in a ratchet model [ij- Thus, HF signals seem to play 
a two-fold role, positive or negative, depending on the 
non-linear system under consideration. 

In this paper we study the influence of a HF field on 
a rather simple excitable model, such as the archetypal 
FHN system. This system is governed by the following 
equations JlS] 

ev = v-v^ -w + ^{t) + S{t)+T{t) (1) 
w = — w ^ h. (2) 

In the context of neurophysiology, v{t) is called the volt- 
age variable, and 'w{t) the recovery variable. Here S{t) 
is an external forcing of period T, ^ {t) a noisy term, and 
T{t) is an added HF signal which will be specified later 



on. The "neuronal" noise £,{t) is assumed to be an unbi- 
ased Gaussian white noise with autocorrelation function 
{£.{t)£,{s)) = 2DS{t — s). The values of the model param- 
eters define the dynamical regimes of the system, and are 
discussed below. 

Let us focus our attention in the deterministic FHN 
model (D = 0) in the absence of HF signal {T{t) = 0). In 
the case of a time-independent external signal S{t) = So, 
the excitable character of the FHN model relies on the 
existence of a threshold value of this signal, ^h, at which 
a Hopf bifurcation takes place. A simple linear stability 
analysis 4] shows that the condition e < 1 is necessary to 
assure the existence of such a threshold Su- In fact, it is 
usually assumed that e <C 1 so that the voltage variable 
is much faster than the recovery variable. In addition, 
we must choose 7 > 1 in order to guarantee that for sub- 
threshold signals there exits a unique stable attractor I3| . 
This is usually the case of interest in most applications of 
the model to real problems, including neuron dynamics. 
In this situation, for Sq < Sn the system is quiescent, 
i.e., the attractor is a globally stable fixed point (■^o, wq). 
When the threshold value Sn is exceeded, the fixed point 
becomes unstable and a stable limit cycle appears. As a 
consequence, when v(t) is viewed as a function of time, 
a sequence of "spikes" is observed, each spike resembling 
an action potential as detected in real neurons. 

In the top-left panel of Fig. ^ we show several rep- 
resentative phase-space trajectories obtained by solving 
the deterministic FHN equations for initial conditions at 
the edges of the frame. The legend p = corresponds 
to r{t) = 0, as will be clarified below. We have used 
the standard parameter values |3j £ = 0.02, 7 = 4, 
and b = 2.8, which corresponds to an excitable regime. 
The external forcing is stationary with = 0.32. Since 
this is a subthreshold value [the threshold value is given 
by Eq. Q with p = 0], all trajectories end up at the 
stable attractor, which is depicted by a solid filled cir- 
cle. With a thick (blue) solid line we have plotted the 
separatrix which separates a region with excitable tra- 
jectories (those with a large excursion in phase space) 
from a region with non-excitable ones (shorter excur- 
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FIG. 1: (Color online) Deterministic time evolution in phase 
space (v, w) for different representative initial conditions and 
a stationary subthreshold external force S(i) — So = 0.32. 
From top-left to bottom-right: p = (no HF signal); p = 0.4; 
p = Pc ^ 0.80829 . . .; and p = 1.2. The stable attrac- 
tors are depicted by filled solid circles. The nuUclines [i.e., 
7W — to -I- 6 = and c{p)v — v^ — w-\-A = 0\ are represented by 
dotted lines. Top panels show the separatrix between nonex- 
citable and excitable behaviors by thick (blue) solid lines. (All 
quantities in dimensionless units.) 

sion) . The excitable paths are reminiscences of the above 
mentioned stable limit cycle, which is only present in 
the superthreshold regime. In the presence of noise, the 
subthreshold dynamics consist of a random trajectory 
around the stable attractor, which eventually crosses the 
separatrix into the excitable region (see Fig. These 
excitations appear as random spikes in the variable w(i), 
as shown in the top-left panel of Fig. 13 Notice that 
in this panel the external forcing S{t) — Acos(r2t) is 
time-dependent. When the external forcing varies slowly 
compared with the characteristic time-scales of the FHN 
model, the above qualitative scheme remains valid if one 
considers a sequence of diagrams corresponding to the 
instantaneous values of S{t). 

Let us now study the effect of a HF signal. In this 
paper, we assume that T{t) takes the form: 

r(i) = NVLr cos{NVtt + ip), (3) 

where £7 = 27r/T, is an arbitrary phase, iV is a large in- 
teger and r is the ratio of the amplitude of the HF signal 
to its frequency. We are interested in situations in which 
the parameters iVJ7r and NVt are much larger than the 
rest of the parameters in the problem. In this sense, we 
will say that T{t) is a strong, high frequency monochro- 
matic force. This situation corresponds formally to tak- 
ing the limit N —* oo while keeping the ratio r fixed. The 
assumption of a strong signal is necessary if the high fre- 
quency signal is to have any effect. 
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FIG. 2: (Color online) Typical phase space trajectory for the 
same parameter values as in the top-left panel of Fig. and a 
noise strength D — 10"''. (All quantities in dimensionless 
units.) 

In this limit it is clear from Eq. ^ that the stochastic 
process v{t) is "fast", since its time derivative is of order 
N . A "slow" stochastic process, w(i), can be defined after 
extracting the "fast" dependence from w(t), 

v{t) = v{t)- p&m{Nnt + ip), (4) 

where p — r/e. Using Eq. in Q)-©, it is easy 
to check that the time derivatives of both v{t) and 
w{t) = w{t) are of order 1 as ^ cxo. Therefore, a 
large number of oscillations of the function psvn.{NQ,t+Lp) 
takes place before a significant change of v{t) and w{t) 
occurs. As a consequence, in this limit, their dynamics 
become independent of the particular value of the phase 
ip, and the resulting equations of motion can be simpli- 
fied by averaging over all possible values of this phase. 
Then, the following equations are obtained 

eb = c{p)v - ^ S{t) ^ S_{t) (5) 
w = ~ w + b. (6) 

where c(p) = 1 - 3pV2, v{t) = (2Tr)-^ J^'" dip v{t), and 
w{t) is defined analogously. 

A linear stability analysis of the deterministic version 
of Eqs. ©-(0) for a time-independent external signal 
^o, shows that the above mentioned condition e < 1, 
which assured the existence of the limit cycle, is re- 
placed in this case by e < c{p). Thus, for all values 
p> Pc = ^2(1 — e)/3 the fixed point is stable regardless 
of the value of So- This fact leads to the disappearance of 
the separatrix for subthreshold signals, as shown in the 
bottom panels of Fig ^ As a consequence, the system 
looses its excitability and no spikes are fired. 

As mentioned before, this description is appropriate 
even in the presence of a slow time-dependent external 
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FIG. 3: Typical trajectories v{t) for a noise strength D — 
5xlO~* and an external force S{t) = Acoa{flt), with A — 0.32 
and n = 0.3. The HF signal is r{t) = Nflpe cos{Nflt), with 
— 4000 and the same values of p as in Fig.0 (All quantities 
in dimensionless units.) 



signal. This case have been considered in Fig. |2| where 
we present four typical trajectories v{t) for the same val- 
ues of p as in Fig. ^ and a slow time-dependent external 
forcing S{t) = Acos{nt), with A = 0.32 and Q = 0.3. 
The remaining parameter values are 13 = 5x10^* and 
N = 4000. The bottom panels show vividly the vanish- 
ment of the spikes due to the suppression of the system 
excitability by the HF signal. 

In addition, the above stability analysis also shows that 
for ratio values p < Pc, the Hopf bifurcation threshold is 
given by 
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Thus, the larger the ratio p, the larger the threshold Sr- 
Consequently, the limit cycle becomes less accessible in 
the presence of the HF signal. This suggests that for sub- 
threshold signals the spike rate decreases as p increases. 
This fact is corroborated by the numerical simulations. 

In order to quantify the effect of the HF signal in the 
FHN model, let us consider the dependence of the spec- 
tral amplification rj on the ratio p. Since we are interested 
in the firing dynamics, first we reduce the stochastic pro- 
cess v{t) to a spike train Q{t) by using the procedure of 
Ref. [13. That is, a rectangular pulse of height unity 
and time length 0.15 is generated whenever v{t) crosses 
a threshold value, chosen as uth — 0.5, having previously 
crossed a reset value = —0.5. We have chosen the 
output variable v{t) given by Eq. instead of v{t), 
to avoid the counting of fake firings due to the rapid 
oscillations. Then, the spectral amplification associated 
to the process 9(t) is defined as the ratio of the inte- 
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FIG. 4: (Color online) Spectral amplification 77 as a function 
of p. The noise strength values are: D = 3x 10"* (diamonds), 
D = 5 X 10"'' (triangles), and D = 7 x 10"* (squares). The 
values of A'^, A, and Q, are the same as in Fig. |3 Solid lines 
depict the results obtained from v {N — > 00). Inset shows rj 
vs. D in the absence of HF signal. (All quantities in dimen- 
sionless units.) 



grated power stored in the first delta-like spike of the 
output power spectral density to that stored in the cor- 
responding delta-like spike of the input power spectral 
density. It can be shown that this quantity is given 
by = A\Mi\^/A^, where Mi = dt {Q{t))e-^^* /T. 

In Fig. ^ we present the dependence of on p for 
three values of the noise strength. Solid lines depict the 
simulation results using Eq. (HJ-® (the case N — > 00), 
whereas points represent the data with N — 4000. No- 
tice that rj vanishes for a given value of p as a conse- 
quence of the disappearance of the model excitability, 
the exact value of p depending on the filter parameters 
vth and Uj-. For the lowest and the intermediate val- 
ues of the noise, £) = 3 x 10^'' and £) = 5 x 10"^, rj 
decreases monotonously with p, although in the second 
case a plateau can be observed for small values of this 
parameter. In contrast, for the largest value of the noise 
strength, D = 7 x 10"*, rj peaks at an optimal value of 
p. This non-monotonic behavior has also been detected 
recently in bistable systems and has been coined vibra- 
tional resonance (VR) [itI H^. However, the situation 
here is just the opposite to that observed in bistable sys- 
tems, because in the deterministic FHN model there is 
no VR effect, and it only appears when the strength of 
the added noise is large enough. 

VR in excitable systems has also been reported in 
Ref. In contrast to our work, in that reference, there 
is no clear-cut separation between the time-scales asso- 
ciated with the HF signal, (7V17)"^, and with the firing 
events. As a consequence, they observe a different be- 
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havior. Besides, in Ref. [l9j, the HF signal is added to 
Eq. (0) instead of to Since u'(t) appears linearly in 
both Eqs. ^ and using the above techniques, it is 
easy to prove that in the formal limit N ^ oo, the HF 
signal would produce no effect. 

In order to understand the behavior observed in Fig. 01 
let us consider first the dependence of rj vs. D in the 
absence of HF signal {p = 0). As shown in the in- 
set of Fig. ^ the spectral amplification displays a non- 
monotonic behavior, with a maximum at an optimal 
noise strength Z3opt at which the coherence of the output 
signal is maximized. This is the celebrated phenomenon 
of stochastic resonance, reported for the FHN model in 
Refs. [ilE3. 

For D < Dopt the coherence of the out- 
put signal is less than the optimal because the number 
of spikes is insufRcicnt to optimize the coherence of the 
output signal. In contrast, for D > Dopt, the number of 
spikes is excessive, and an appreciable number of them 
are fired out of phase with respect to S{t), which de- 
grades the coherence. Diamonds in Fig. 0] correspond to 
a noise strength less than the optimal. Since the pres- 
ence of the HF signal decreases further the number of 
spikes, rj must decrease monotonously as p increases. On 
the other hand, when D > Z?opt (squares in Fig.^J), the 
decrease in the firing rate produced by the HF signal re- 
moves the incoherent extra spikes. This fact brings the 
system near the optimal operation regime, which results 
in a higher value of rj than the corresponding to p = 0. Of 
course, when p is increased further, the number of spikes 
decrease until total disappearance at a critical value of 
p, as discussed above. Consequently, a non-monotonic 
behavior of as a function of p is expected for a fixed 
value of _D > Dopt- Triangles in Fig. ^correspond to an 
intermediate value of D, close to the optimal. 

In conclusion, we have shown that while a HF signal is 
able to improve the response of the FHN neuron system 
when the noise strength is larger than Dopt, when p is 
large enough, the HF signal inhibits completely the firing 
activity. This suggests a negative effect in the neuron 
activity for sufficiently strong HF signals which could be 
related to the functional biological damages caused by 
electro-pollution \vi\ . 
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